Multiple Plane Phase Retrieval Based On Inverse Regularized Imaging and Discrete Diffraction Transform

The phase retrieval is formulated as an inverse problem, where the forward propagation is defined by Discrete Diffraction Transform (DDT) [1], [2]. This propagation model is precise and aliasing free for pixelwise invariant (pixelated) wave field distributions in the sensor and object planes. Because of finite size of sensors DDT can be illconditioned and the regularization is an important component of the inverse. The proposed algorithm is designed for multiple plane observations and can be treated as a generalization of the Gerchberg-Saxton iterative algorithm. The proposed algorithm is studied by numerical experiments produced for phase and amplitude modulated object distributions. Comparison versus the conventional forward propagation models such as the angular spectrum decomposition and the convolutional model used in the algorithm of the same structure shows a clear advantage of DDT enabling better accuracy and better imaging.